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Abstract. In the 2-dimensional n-body problem, n > 3, in spaces of 
constant curvature, k 0, we study polygonal homographic solutions. 
We hrst provide necessary and sufficient conditions for the existence 
of these orbits and then consider the case of regular polygons. We 
further use this criterion to show that, for any n > 3, the regular n- 
gon is a polygonal homographic orbit if and only if all masses are equal. 
Then we prove the existence of relative equilibria of non-equal masses 
on the sphere of curvature k > for n = 3 in the case of scalene 
triangles. Such triangular relative equilibria occur only along fixed 
geodesies and are generated from fixed points of the sphere. Finally, 
through a classification of the isosceles case, we prove that not any 
three masses can form a triangular relative equilibrium. 

1. Introduction 

We study here the curved n-body problem, defined as the motion of n 
point particles of masses mi, rri2, ■ ■ ■ , m n > in spaces of constant curva- 
ture, k^0, under the influence of a natural extension of Newton's gravi- 
tational law. The potential defining this law is given by the force function 
U K , whose expression appears in (J3J). The corresponding potential, given 
by — U K , preserves the basic properties of its Euclidean analogue: it is 
a harmonic function in the 3-dimensional space, i.e. it satisfies Laplace's 
equation, and generates a central field in which bounded orbits are closed, 
[21] , in agreement with an old result proved by Joseph Louis Bertrand for 
the Euclidean case, [2]. 

The research direction we are following here started in the 1830s, when 
Janos Bolyai and Nikolai Lobachevsky independently proposed a curved 
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2-body problem in the hyperbolic space H 3 , given by a force proportional 
to the inverse of the area of the sphere of radius equal to the distance 
between bodies. This problem was studied by top mathematicians, such 
as Lejeune Dirichlet, Ernest Schering, [31], [32], Wilhem Killing, [18], |19J . 
[20] . and Heinrich Liebmann, [23], [21], [25]. Schering was the one who 
came up with an analytic form of the potential, which can be expressed 
in terms of the cotangent of the distance, for k > 0, and the hyperbolic 
cotangent of the distance, for k < 0. The Newtonian law is recovered 
when k — y 0. Recently, Jose Carinena, Manuel Ranada, and Mariano 
Santander, [3J, proved several new results for this 2-body problem and 
revisited some old properties with the help of modern methods. The study 
of the quantum analogue of the curved 2-body problem was proposed by 
Erwin Schrodinger, [33] . and continued by Leopold Infeld, [IB] , and Alfred 
Schild, [TZ]. 

Other attempts at extending the Newtonian case to spaces of constant 
curvature, such as the efforts of Rudolph Lipschitz, [26], did not sur- 
vive, mostly because the proposed potentials lacked the basic physical 
properties mentioned above. Unlike Liebmann, who showed that all the 
fundamental orbits of the Kepler problem have analogues in curved space, 
Lipschitz could not obtain explicit solutions; he only succeeded to express 
the orbits in terms of elliptic integrals. 

In the direction of research proposed by Bolyai and Lobachevsky, the 
generalization of the equations of motion to any n > 2 was given in [Z] , a 
paper posted in arXiv, but submitted for publication as two connected pa- 
pers, [8] and [9], in which we obtained a unified framework of approaching 
the problem for any k^O. We also proved there the existence of several 
classes of relative equilibria, including the Lagrangian orbits, i.e. the case 
n = 3 of the equilateral triangle. Relative equilibria are orbits for which 
the configuration of the system remains congruent with itself for all time, 
so the mutual distances between bodies are constant during the motion. 

It is worth mentioning that the study of the curved n-body problem, 
for n > 3, might help us better understand the nature of the physical 
space. Gauss allegedly tried to determine the geometry of the universe by 
measuring the angles of a triangle formed by the peaks of three mountains. 
Even if the goal of his topographic measurements was different from what 
anecdotical history attributes to him (see [25]), this method of deciding 
the nature of space remains valid for astronomical distances. But since 
we cannot measure the angles of cosmic triangles, we could alternatively 
check whether certain (potentially observable) celestial motions occur in 
the universe, and thus decide whether the physical space has negative, 
zero, or positive curvature. 
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Specifically, we showed in [7J and [8] that while Lagrangian orbits of 
non-equal masses are known to occur for k — 0, they must have equal 
masses for n ^ 0. Since Lagrangian solutions of non-equal masses exist 
in our solar system (such as the triangle formed by the Sun, Jupiter, and 
the Trojan asteroids), we can conclude that, if assumed to have constant 
curvature, the physical space is Euclidean for distances of the order 10 1 
AU. The discovery of new orbits of the curved n-body problem might help 
us better understand the large-scale geometry of the universe. 

The most recent papers on the curved n-body problem deal either with 
singularities, [5], [DJ, a subject we will not approach here, or with homo- 
graphic solutions and, in particular, with homothetic orbits and relative 
equilibria, [7J, [9j, [6j. Homographic solutions are orbits whose configu- 
ration remains similar to itself all along the motion. In particular, when 
rotation takes place without expansion or contraction, the homographic 
orbits are called relative equilibria, as mentioned earlier. They behave 
like rigid bodies, maintaining constant mutual distances. Homothetic so- 
lutions are homographic orbits that experience expansion and/or contrac- 
tion, but no rotation. 

The homograpic solutions can be put in a broader perspective. They are 
also the object of Saari's conjecture, [3D], [TDJ, which we partially solved 
for the curved n-body problem, [7J, [Sj, as well as Saari's homographic 
conjecture, [30], [11]. Both have recently generated a lot of interest in 
classical celestial mechanics (see the references in [TDJ, [TTJ) and are still 
unsolved in general. 

In the classical Newtonian case, [M], as well as in more general clas- 
sical contexts, [I], the standard concept for understanding homographic 
solutions is that of central configuration. We will not employ it here since 
most computations appear to be simpler without using it. The reason 
for these complications is connected to the absence of the integrals of 
the centre of mass and linear momentum from the curved n-body prob- 
lem. These integrals seem to be specific only to Euclidean space. Indeed, 
n-body problems derived by discretizing Einstein's field equations, as ob- 
tained by Tullio Levi-Civita, [21] , [22], Albert Einstein, Leopold Infeld, 
and Banesh Hoffmann, [T2], and Vladimir Fock, [T3], also lack such inte- 
grals. 

In this paper we study polygonal homographic orbits of the 2-dimension- 
al curved n-body problem. In Section 2, we introduce the notation and 
the equations of motion as well as their first integrals. In Section 3 we 
define polygonal homographic orbits and their basic particular cases: ho- 
mothetic orbits and relative equilibria. We also provide a motivation for 
our definitions, which are given in the spirit of the Euclidean case. In 
Section 4, we state and prove necessary and sufficient conditions for the 
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existence of polygonal nomographic orbits. Section 5 is dedicated to the 
study of regular polygons. We show, on one hand, that if the masses are 
equal, the regular n-gon is a solution of the equations of motion for any 
n > 3. On the other hand, we prove that regular n-gons can be solutions 
only when the masses are equal. These results extend the Perko-Walter- 
Elmabsout theorem, [29], [13], to spaces of non-zero constant curvature. 
In Section 6, we focus on the case n = 3 and prove that the equations 
of motion admit no homographic orbits if the triangle is not equilateral. 
For k > this result is true as long as the homographic orbit is not a 
relative equilibrium rotating along a great circle of the sphere. In Sec- 
tion 7 we discuss the orbits omitted in the previous section, namely the 
relative equilibria that move along a geodesic of the sphere. We prove 
the existence of relative equilibria of non-equal masses in the case of sca- 
lene triangles. Such triangular relative equilibria occur only along fixed 
geodesies and are generated from fixed points on the sphere. Finally we 
show that not any three masses can form a triangular relative equilibrium 
by providing a large class of counterexamples in the isosceles case. 



In this section we introduce the equations of motion of the curved n- 
body problem on 2-dimensional manifolds of constant curvature, namely 
spheres embedded in IR 3 , for n > 0, and the upper sheets of hyperboloids 
of two sheetg^ embedded in the Minkovski space IR 2,1 , for n < 0. 

Consider the point particles (bodies) of masses mi, m2, ■ ■ ■ , m n > 
in IR 3 , for n > 0, and in IR 2,1 , for k < 0, whose positions are given 
by the vectors q { = (x iy zi), i = 1, n. Let q = (qi, q 2 , . . . , q„) be the 
configuration of the system, and p = (pi, p 2 , . . . , p n ), with p,; = m^, % = 
representing the momentum. We define the gradient operator with 
respect to the vector q» as 



2. Equations of motion 




where a is the signum function, 
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for k > 
for k < 0, 




,,Vq n ). For the 3-dimensional vectors a 
b z ), we define the inner product 



The upper sheet of the hyperboloid of two sheets corresponds to Weierstrass's 
model of hyperbolic geometry (see Appendix in [7] or [5]). 
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and the cross product 

(3) a <g> b := [a y b z - a z b y , a z b x - a x b z ,a(a x b y - a y b x )). 

The Hamiltonian function of the system describing the motion of the 
n-body problem in spaces of constant curvature is 

H K (q, p) = T K (q, p) - U K (q), 

where 

1 n 

r«(q,p) = ^ ^2 m 7\Pi ©Pi)(«qi 0q*) 
i=i 

defines the kinetic energy and 

(A) II (a) V m^l^l^/tq.Qq, 

l<fen M** ® ° ^ " ^ ° ^) 2 ] 1/2 

is the force function, — U K representing the potential energjU. Then the 
Hamiltonian form of the equations of motion is given by the system 



(5) 



Pi = V qi C/«(q) - m i 1 n{p i pi)qi, i = 1, n, k ^ 0, 
where the gradient of the force function has the expression 

m v /7 r ^ = V m « m jl K l 3/2 ( K qj Q qj)[(^ Q q»)qj ~ ( K q» Q gjjqj 

u 1 ^ K^Oq,)^©^-)-^©^) 2 ] 372 ' 

The motion is confined to the surface of nonzero constant curvature k, 
i.e. (q,p) e T*(M 2 ) n , where T*(M 2 ) n is the cotangent bundle of the 
configuration space (M 2 ) n , and 

M 2 K = {(x, y, z)el 3 k(x 2 + y 2 + az 2 ) = 1}. 

In particular, M 2 = S 2 is the 2-dimensional sphere, and M 2 ^ = H 2 is 
the 2-dimensional hyperbolic plane, represented by the upper sheet of the 
hyperboloid of two sheets (see the Appendix of [7j or [8] for more details). 
We will also denote M 2 by S 2 for k > and by H 2 for k < 0. 

Notice that the n constraints given by ncy, © q« = 1, i = 1, n, imply that 
q« © P« = 0, so the 6n-dimensional system (jSJ) has 2n constraints. The 
Hamiltonian function provides the integral of energy, 

#«(q, P) = h, 



2 In [7] and [3], we showed how this expression of U K follows from the cotangent 
potential for k ^ 0, and that Uq is the Newtonian potential of the Euclidean problem, 
obtained as n — > 0. 
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where h is the energy constant. Equations (EJ) also have the three integrals 
of the angular momentum, 



n 



(7) 2^ Qi ® Pi = c, 

i=l 

where c = (a,/?, 7) is a constant vector. Unlike in the Euclidean case, 
there are no integrals of the center of mass and linear momentum. Their 
absence complicates the study of the problem since many of the standard 
methods don't apply anymore. 

Using the fact that /tq^ q, = 1 for i — l,n, we can write system ([5]) 

as 



qi = > -^7 r t^& 2 («qi © qjjqi, i = 1, n, 

[a - a(KqiQ q j ) 2 \ 



which is the form of the equations of motion we will use in this paper. 
The sums on the right hand side of the above equations represent the gra- 
dient of the potential. When k — > 0, both the sphere (for k — > 0, n > 0) 
and the upper sheet of the hyperboloid of two sheets (for k — > 0, k < 0) 
become planes at infinity, relative to the centre of the frame. The seg- 
ments through the origin of the frame whose angle measures the distance 
between two points on the curved surface become parallel and infinite, so 
the distance in the limit plane is the Euclidean distance. Consequently 
the potential tends to the Newtonian potential as k — > (see [7] or [8] 
for more details). The terms involving the velocities occur because of the 
constraints imposed by the curvature. They vanish when k — > 0. 

3. Polygonal homographic orbits 

In this section we define the polygonal homographic solutions of the 
curved n-body problem as well as two remarkable subclasses of solutions: 
the polygonal homothetic orbits and the polygonal relative equilibria. 
Then we justify the content of these definitions. 

Definition 1. A solution of equations (151), which describe the curved 
n-body problem, is called polygonal homographic if the bodies of masses 
m 1; m 2 , . . . , m n > 0, with n > 3, form a polygon that is orthogonal to the 
z axis and remains similar to itself for all time t. 

According to Definition [H the size of a polygonal homographic solution 
can vary, but its shape remains the same. Notice that we imposed the 
condition that the plane of the polygon is always perpendicular to the z 
axis. This condition is equivalent to saying that all masses have the same 
coordinate z(t), which may vary in time. This condition is not imposed 
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for mere simplicity but because polygonal nomographic solutions may not 
exist without it. Though a complete proof of this conjecture is still eluding 
us, we will explain later in this section why we think this property is true. 

We can represent a polygonal homographic solution of the curved n- 
body problem in the form 

(9) q = (qi, • • • ,q„), q* = (x u y i: Zi), 

Xi = r cos(u + at), y% = r sin(w + at), z% = z, i = 1, n, 

where < ot,\ < < ■ ■ ■ ct n < 2ir are constants; the function z = z(t) 
satisfies z 2 = an" 1 —err 2 ; a is the signum function defined in (JTJ; r := r(t) 
is the size function; and to := u(t) is the angular function. 

Indeed, for every time t, we have that xf{t) + yf(t) + crzf(t) = i = 
1, n, which means that the bodies move on the surface M 2 ., and the angles 
between any two bodies, seen from the centre of the circle containing 
the polygon, are invariant in time. Therefore representation fl9]) of the 
polygonal homographic orbits agrees with Definition [TJ 

Definition 2. A polygonal homographic solution of equations (jHJ), which 
describe the curved n-body problem, is called polygonal homothetic if the 
polygon having at its vertices the bodies of masses mi,m 2 , . . . ,m n > 0, 
with n > 3, expands or contracts, but does not rotate around the z axis. 

In terms of representation 0, a polygonal homographic solution is 
polygonal homothetic if u(t) is constant, but r(t) is not. Such orbits occur, 
for instance, when n bodies of equal masses, lying initially at the vertices 
of a regular polygon inscribed in a non-geodesic circle of the sphere, are 
released with zero initial velocities, to end up in a total collision. 

Definition 3. A polygonal homographic solution of equations ([8]) ; which 
describe the curved n-body problem, is called a polygonal relative equilib- 
rium if the polygon having at its vertices the bodies of masses mi, ni2, . . . , 
m n > 0, with n > 3, rotates around the z axis and maintains fixed mutual 
distances. 

In terms of representation a polygonal relative equilibrium occurs 
when r(t) is constant, but u(t) is not. These orbits have a rich recent his- 
tory and have been extensively studied in the general context of geometric 
mechanics (see, e.g., [I], [27]). 

We will further loosely use the terms "dynamical polygon" or "dynam- 
ical n-gon" to describe any of the polygonal homographic, homothetic, or 
relative equilibrium orbits that will occur in this paper. These terms will 
also occur when we check potential solutions, prior to knowing whether 
they satisfy the equations of motion. 
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Using the concept of relative equilibrium, we can now provide a justifi- 
cation for introducing the orthogonality condition in Definitions 1, 2, and 
3. Consider a polygon formed by the bodies of masses m 1; m 2 , . . . , m n > 0, 
which move on the surface Mj* according to equations (J5J) . Unlike in repre- 
sentation (Q, the plane of the polygon can have any angle (not only 7r/2) 
relative to the z axis. We assume this angle to be constant in time. We as- 
sign initial velocities such that the dynamical polygon rotates around the 
z axis, which passes at all times through the same point inside the poly- 
gon. The Principal Axis Theorem (see Appendix in [7J or [S]) guarantees 
the validity of this scenario without any loss of generality. 

Let us now seek necessary conditions for the existence of the above de- 
scribed motion. Notice that the projection of the polygon to the xy plane 
is, at any time t, a polygon congruent with the projection obtained at 
t — 0. We can then describe the motion of the real polygon in terms of the 
angles of the projected polygon. So let us assume that the relative equilib- 
rium is represented by the coordinates qi(t) = (xi(t), yi(t), Zi(t)), i = l,n, 
given by 

(10) Xi(t) = Ti cos(fit + ttj), yi(t) = Ti sin (fit + otj), Zi(t) = Z{ (constant) 

where < a.\ < a 2 < ■ ■ ■ a n < 2n, fi ^ 0, and r 4 > are constants; 
zf = ctk -1 — arf; and a is the signum function defined in (Tj[|). 

For such a solution to exist, it is necessary that the total angular mo- 
mentum is the zero vector or a vector parallel with the z axis. Otherwise 
the angular-momentum vector would rotate around the z axis, in viola- 
tion of the angular-momentum integrals (J7j). This means that at least the 
first two components of the vector Y^=i m i ( \i ® q« are zero, i.e. 

n n 

^rrii^yiZi - ijiZi) = y~]mi(xjZi - XiZi) = 0. 

i=l i=l 

Using ffTQ]) and the fact that fi ^ 0, these two equations take the form 

n n 

raiTiZi cos(fit + «j) = raiTiZi sin(fit + ai) = 0. 
i=i i=i 

In general, the above conditions for the existence of motions described 
above are not satisfied for all values of t. But there are exceptions, such 
as when all Z{ (and implicitly all Tj) are equal, all m ; are equal, and 
the polygon is regular. So the case when all z% are equal, which implies 
orthogonality relative to the z axis, seems like a good point to start from. 
This position is supported by the proof that Lagrangian solutions must be 
orthogonal to the z axis (see [7] or [8]). These remarks justify our choice 
of the orthogonality condition in the definitions of this section. 
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4. Necessary and sufficient conditions 

The goal of this section is to state and prove two equivalent criteria that 
provide necessary and sufficient conditions for the existence of polygonal 
homographic solutions of the curved n-body problem. The first criterion 
can be expressed as follows. 

Criterion 1. Consider n > 3 bodies of masses m 1; m 2 , . . . , m n > mov- 
ing on the surface M 2 .. The necessary and sufficient conditions that a 
solution of the form is a polygonal homographic orbit of equations OH]) 
are given by the equations 

(11) S 1 = S 2 = ... = S n and 7l = 72 = • • • = 7n , 

where 

n n 

(12) Sj = } j rrijfiji, 7i = ^ rrijVji, i = l,ra, 



V 7 Pjl 1/2 Co 2^/2' ^ 3 / 2 Co 2^/2' 

Cj-j (2 - c^nr*) 6 ' 2 Cj- (2 - CjiKr z ) 6 ' A 

(14) = sin (a,- - cji = 1 - cos(aj - c^), z, j = 1, n, i ^ j. 



Proof. Let us check in equations (JS]) a solution of the form ([9]). For this 
purpose we first compute that 

ii = r cos(o; + ai) — ruj sin(w + ai), iji = r sin(o; + ai) + ru cos(u + ai) 

Zi = —arr(aKT 1 — ar 2 )^ 1 ^ 2 . 
Xi = (r — rcu 2 ) cos(u; + ai) — (riu + 2fu) sm(co + ai), 
yi = (r — ru 2 ) sm(u + ai) + (ru + 2ru) cos(u + ai), 
'ii = -arr\aK~ x - ar 2 )~ l/2 - n~ l r 2 (aKT l - ar 2 )~ 3/2 . 

Some long but straightforward computations for the Zi component lead 
us to to the equations 

. o 

(15) f = r(l - nr 2 )u - - - A i; i = 1, n, 

1 — KT Z 

where 

, , \ mj 1 — fix 

(!6) **= £ "T/^ b? * = 1 ^ 

and the constants Cjj are defined in (fT4l) . Notice that, as long as the bodies 
stay away from collisions, we have Cji > 0. 
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From the equations corresponding to and jji, we further apply the 
formula for the cosine of the sum of two angles to cos(u; + atj), viewed as 
cos[(w + cxi) + (atj — «j)]. Then, using the fact that the equations must 
be satisfied for all times t, and comparing the similar terms, we recover 
equations (TT5T) as well as obtain some new equations, 

(17) ruj + 2ruj — Tj = 0, i = l,n, 

where 

n 

(is) r,= 3/2 — r> * = T ^ 

and the constants Sji,Cji, i,j = l,n,i ^ j, are defined in (fill) . Notice 
that Sji can have any sign, with Sji = only if a 3 - — a.{ = ±7r/2 or ±37r/2. 

Equations ( 1T51) and ( 1TTI) describe the motion of the polygonal homo- 
graphic orbit with the help of the size function, r, and the angular func- 
tion, to. But this system of differential equations makes sense only if the 
conditions 

Ai = A 2 = • • ■ = A„ and T 1 = T 2 = ■ ■ ■ = T n 

are satisfied, where the functions A^r^i = l,n, are defined in (|T6j) and 
(|T8|) . respectively Then, for any initial conditions (r(0), w(0)) that are not 
singular (i.e. are not collisions for any k^O, and are not antipodal for k > 
0, see j?j or [S]), we are guaranteed the local existence and uniqueness of 
an analytic solution for the system given by equations (fT5|) and (fT7|) . Since 
the phase space is a connected manifold, this solution can be analytically 
extended to a global solution, defined either for all time or until the orbit 
reaches a singular configuration, as it may happen, for instance, in the 
case of a homothetic orbit that ends in a total collision. 

Multiplying equation (fl5j) by r 2 /(l — nr 2 ) and equation (fT7|) by r 2 , the 
conditions for the existence of polygonal homographic orbits reduce to 

S 1 = 5 2 = ■ ■ ■ = 5 n and j ± = j 2 = ■ ■ • = 7„, 

as defined in the above statement. This remark completes the proof. □ 

Remark. Notice that, if S denotes any 5i and 7 denotes any 7$, % = l,n, the 
functions fiji and 6 are always positive, while Vji and 7 can be negative, 
positive, or zero. 

We can restate Criterion [TJ in terms of linear algebraic systems as fol- 
lows. The equivalence between Criterion [T] and Criterion |5] is obvious. 

Criterion 2. Consider n > 3 bodies of masses m 1; m 2 , . . . , m n > mov- 
ing on the surface M 2 .. The necessary and sufficient conditions that a 



Polygonal nomographic orbits of the curved n-body problem 



11 



solution of the form (jHJ) is a polygonal homographic orbit of equations (jHJ) 
are that there exist 5 > and 7GI such that the linear systems 

(19) Am T = d T and rm T = e r 

have the same set of solutions, where 

m = (m 1 ,m 2 ,...,m n ), d = ( 5,6, .,5) , e = ( QVL - -_rf) , 

n times n times 



A 



/!21 A*31 

A*i2 /i 32 

A*13 A*23 

Mln /^2n /^3n 



/^n2 
A*n3 





V 2 \ vz\ 

V\2 Z/32 
^13 ^23 







. V\ n U2n ^3n 







i ^ j are defined in ([13]) , and the upper index T denotes 
the transpose of a vector or matrix. 

Remark. Notice that A is symmetric, i.e. A T = A, whereas T is skew- 
symmetric, i.e. T T = — r. In the next section we will see that both A and 
T are circulant matrices, i.e. each row vector is rotated one element to the 
right relative to the previous (upper) row vector. 



5. Regular polygons 

In this section we will study the case of regular polygons and prove 
that the curved n-body problem, with n > 3, admits regular polygonal 
homographic orbits if and only if all the masses are equal. This result 
is valid only for n > 4 in the Euclidean case, where it is known as the 
Perko-Walter-Elmabsout theorem. The classical Lagrangian orbits, given 
by the dynamic equilateral triangle, allow non-equal masses, a situation 
that seems to be specific to the Euclidean space, most likely because 
some symmetries are lost if the curvature is not zero (see also jl]). To 
emphasize certain details of our result, we will separately consider the 
direct and converse components of this theorem. 

Theorem 1. Consider the curved n-body problem, n > 3, given by system 
OH]). If n bodies of equal masses, m := mi = vrt2 = ■ ■ ■ = m n , lie initially 
at the vertices of a regular n-gon parallel with the xy plane, then there 
is a class of initial velocities for which the corresponding solutions are 
homographic. These orbits also satisfy the equalities 71 = 72 = • • ■ = 

7n = 0. 

Proof. Let the dynamical n-gon be regular and, for every time instant, lie 
in a plane parallel with the xy plane. We will show that 8% = 82 = ■ • • = S n 
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and 71 = 72 = • • • = 7n, identities which, by Criterion CEJ assure the 
existence of the desired homographic solution. 

Without loss of generality, we assume that a iy i = l,n, the constants 
which provide the representation (Q of the homographic orbits, are 

cci = 0, «2 = 2ii/n, «3 = 47r/n, . . . , a n = 2(n — l)ir/n. 



Then the differences ay — ai, i,j = l,n, with i ^ j, can take only the 
values 

±2n/n, ±4vr/n, . . . , ±2(n - l)n/n. 

Consequently, these angles fully determine the constants Sij and Cij defined 
in (03). 

We will first study the matrix A. Let s = [n/2] denote the integer part 
of n/2. Then, for n odd, the matrix A takes the circulant form 






ai 


0,2 


a.s-1 


a s 


«s 


a s -i 


a 3 


(12 


ai 


"1 





ai 


a.s-2 


a s -i 


«s 


a s 


(24 


az 


(1-2 


0,2 


01 





a s -3 


a s -2 


a s -i 


a s 


as 


('4 


0-3 


ds-1 


a-s-2 


a s -3 





ai 


(12 


a 3 


a s _i 


(Is 


(Is 


a s 


a.s-1 


a.s-2 


ai 





ai 


a.2 


a s -2 


a a -l 


(is 


a s 


a s 


a.s-1 


a 2 


"i 





ai 


a-s-3 


a-s-2 


a a -i 


a s -i 


«s 


a s 


a 3 


(12 


ai 





a 3 _4 


a s -3 


a s -2 


a 3 


0,4 


a 5 


a.s-1 


a s -2 


a s -3 


a s _4 





ai 


0,2 


(12 


o:i 


«4 


a B 


a s -i 




a 3 -3 


ai 





(11 


«1 


a 2 


a 3 


Cls 


(Is 


a s -i 


a s -2 


a 2 


«i 






whereas for n even, A becomes the circulant matrix 



o 

a2 ai 



ai 




((_> 
fi 




a s a s 
i s -i a 
is-2 a B 



«4 



(12 

a-3 

(14 



ai 

«2 
(1 3 



-2 a s _3 
-l a s _2 
j a s _i 




ai 

(12 



ai 

ai 



«2 

ai 




« 3 
(12 
(II 



(14 
'-'3 
"2 



'-'5 
(14 
(13 



a s -i a s _2 




ai 

(12 



«1 



"1 



(12 
(II 





where ai, . . . , a s > represent the fiji functions in their corresponding 
positions. Notice that in the ith row, the functions fiji with \i— j\ = s and 
\z — j\ — n — s, s — 1, [n/2], are equal, and they are independent of i. But 
in every row of A, there are exactly two functions fiji with \z — j\ — s or 
\i — j\ = n — s, s = 1, [n/2]. For n even, the single /ijj term not captured 
in the above description corresponds to aj — a,, = ir and is of the form 



1/[4(1 



nr 



! ) 3 ' 2 ], so it is also independent of i. Consequently each 
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row of A contains the same elements, only differently ordered. Since all 
the masses are equal, we can conclude that Si = 5 2 = ■ ■ ■ = S n . 
For n odd, the matrix T is circulant, 






6l 


6 2 


bs-i 


bs 


-b s 


bs-i 




-63 


-62 




-bi 


-bi 





6l 


b a -2 


b s -l 


b s 


-b s 




-b 4 


-b 3 




-62 


-6 2 


-61 





b B -3 


b s -2 


b B -i 


b a 




-6s 


-bi 




-63 




b s -l 


-b a -2 - 


bs-3 ■■■ 





6l 


l>2 


63 




bs— 1 


b a 




-6s 


-b s 


'bs-i - 


-b a -2 ■■■ 


-61 





bi 


62 




b s — 2 


bs—i 




% 


b s 


-b a 


-bs-i ■■■ 


-62 


-bl 





61 




bs-3 


b s -2 


bs-i 


ba-l 


b s 


~b 3 


-63 


-b 2 


-61 









f) s g 


bs-2 


&3 


In 


65 


— b s -l 


-bs-2 


-bs-3 


bs-4 







bl 


62 


l>2 


bs 


64 


-6 S - 


-b B -i - 


-b s -2 - 


6 S _3 




-bi 





bi 


6l 


62 


63 




-b s 


-bs-i 


bs-2 




-62 


-bi 







lereas for n even, T takes the circulant form 



















'.'1 


b 2 


6„_i 





-b a -i 




-63 


-62 


-bi 








-61 





h 


b s -2 


b s -i 







-64 


-63 


-62 








-62 


-61 





b s -3 


b s -2 


b a -i 




-bs 


-64 


-63 








-&s-l 


-b s -2 


-b s -3 





bi 


62 




bs-2 


bs-i 














-b s _i 


-b s -2 


-h 





'.'1 




bs-3 


bs-2 


bs-i 




? 




b s -i 





-b 3 -i 


-62 


-61 







bs-4 


bs-3 


bs-2 








b 3 


&4 


bs 


-b s -2 


-bs-3 


-b s -4 







bi 


62 








b 2 


63 


64 


-6 a _i 


-b B -2 


-b B -3 




-bi 





bi 








. &i 


62 


b 3 





-bs-i 


-bs-2 




-62 


-bi 










where 61, 62, • • • , b s > represent the corresponding uji functions. Thus, 
for T, we have to slightly amend the arguments used for A: the two 
functions Vji with \i — j\ = s or \i — j\ = n — s, which are independent 
of i and occur in every row, are equal only in absolute value; they have 
opposite signs because Sji := sin(<x,-— c^) is odd. For n even, the extra term 
corresponding to otj — ai = re is b s = 0. Therefore 71 = 72 = • • ■ = 7« = 0. 
By Criterion [TJ we can now conclude that the dynamical regular ra-gon 
having equal masses at its vertices is a homographic orbit of the curved 
n-body problem. This remark completes the proof. □ 

We can now state and prove the converse of Theorem [TJ 

Theorem 2. // the masses mi, . . . ,m n > 0, n > 3, form a polygonal 
homographic solution of the curved n-body problem given by equations 
such that the polygon is regular, then nil = m 2 = ■ ■ ■ = m n . 

Proof. Notice first that, from the form of the matrix T represented above 
(for both n even and odd), if a dynamical regular n-gon parallel with the 
xy plane is a homographic solution of the curved n-body problem, then 
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the masses must satisfy the system rm T = T . By Criterion [2], the masses 
must also satisfy the system Am T = d T . 

Since, in general, the functions and Uji vary in time, and the solu- 
tions mi, m 2 , . . . , m n of the systems Am T = d T and rm T = T depend 
on and z^j, only solutions for which the masses are constant and posi- 
tive lead to homographic solutions. Therefore if we fix an arbitrary time t 
and show that the corresponding systems with constant coefficients have 
solutions only when vii\ = m 2 = ■ ■ ■ = m n , then no other solutions are 
possible for those systems when t varies. By Theorem [H we can then 
conclude that mi = m 2 = ■ ■ • = m n is the only case when the systems 
with variable coefficients have solutions. 

So let us assume t fixed and start with the system Am T = d T . We 
already know from Theorem [T] that this system has infinitely many solu- 
tions, namely mi = m 2 = • • • = m n = a for any a > 0, where a is some 
function of 5 > 0. But if we fix a value of 5, then a is also fixed. 

We will now show that det A ^ 0. To prove this fact, notice that, except 
for the diagonal elements, all the elements of A are positive. Therefore 
A is positive definite. Indeed, a simple computation shows that for any 
nonzero vector z = (zi, z 2 , ■ ■ ■ , z n ) G M n , zAz T > 0. By Sylvester's crite- 
rion (see, e.g., |T3] ) , positive-definite matrices have positive determinants. 

According to Cramer's rule, the linear system Am T = d T has a unique 
solution, which must be none else than m\ = m 2 = ■ ■ ■ = m n = a, 
with a > fixed. By Theorem [TJ this solution also satisfies the system 
rm T = T . Since the choice of a depends on the choice of 5, and 5 > 
is arbitrarily fixed, we can draw the same conclusion for any choice of 5. 
Consequently a dynamical regular n-gon can be a homographic solution 
of the curved n-body problem only if all masses are equal. This remark 
completes the proof. □ 



6. Non-geodesic scalene triangles 

It is natural to ask whether irregular polygons could form homographic 
orbits. The conditions of Criterion [T] for the existence and uniqueness of 
dynamical n-gons suggest that, in general, this is not the case. Indeed, 
system flTTT) has 2n — 2 linear equations and n unknowns: m 1 , m 2 , . . . , m n . 
Even the case n = 3 leads to a linear system of 4 equations with 3 un- 
knowns, which is still unlikely to have solutions in general. 

We will next prove that homographic orbits of the curved 3-body prob- 
lem cannot exist for | z \ 7^ if the triangle is not equilateral. In other 
words, the symmetries of the equilateral triangle (and consequently the 
equality of the masses) are a necessary condition for the existence of such 
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orbits. The restriction \z\ 7^ is necessary only for k > (being automati- 
cally satisfied for k < 0) because equations (fT5|) do not allow homographic 
orbits to pass through the equator of the sphere, since the function 1 — nr 2 , 
which appears as a denominator, cancels in that case. However, as we will 
see in the next section, non-equilateral relative equilibria moving along the 
equator do exist. 

We can now state and prove the following result. 

Theorem 3. Consider the curved 3-body problem, given by equations (JS} 
with n = 3 and masses m 1 ,m 2 ,m 3 > 0. These equations admit no ho- 
mographic orbits given by scalene non- equilateral triangles for k < 0. For 
K > , they don't admit such solutions either if the bodies stay away from 
the equator z — 0. 

Proof. Notice first that any homographic orbit of the 3-body problem 
must correspond to an acute triangle. Otherwise, at every time instant, 
there is a plane containing the z axis such that all three bodies are on 
one side of the plane (at most two of them in the plane). Therefore the 
total angular momentum vector at that time instant cannot be zero or 
parallel with the z axis. But after rotating by it radians, the bodies reach 
a position on the other side of the plane, and the angular momentum 
has certainly a different direction, in violation of the angular-momentum 
integrals ([7]). 

For the masses m 1 ,m2,m 3 > 0, the conditions in Criterion [2] can be 
written as 

(20) Am T = d T and rm T = e T , 

m = (mi,m 2 ,m 3 ), d = (5,5,5), e= (7,7,7), 





" 


a 


b ' 







u 


V 


A = 


a 





c 


, r = 


—u 





w 




b 


c 







—v 


— w 






a := fi 2 i = A*i2, b := /i 3 i = yUi 3 , c := /i 32 = /i 2 3, 

u := ^21 = -V12, v := u 3 i = -u 13 , w := z/ 32 = -u 23 , 

and fiji, Uji, i,j — 1, 2, 3, % ^ j, as defined in ffTBT) . 

Multiplying the first equation of the system rm T = e T by —w, the 
second by v, and adding them, we obtain the equation 

—uvnii — uwm,2 = 7(1* — w). 

Multiplying the third equation of that system by — u and adding the 
above equation to it, we get the condition j(v — u — w) = 0. Circular 
permutations lead us to the system of conditions 

j(v — u — w) = 7(10 — v — u) — j(u — w — v ) = 0. 
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This system is satisfied either if 7 = or when 

v — u — w = w — v — u = u — w — V — 0. 

But the above equations imply that u — v — w — 0, which is a solution 
with no dynamical consequences, so necessarily 7 = 0. Consequently the 
linear system rm T = e T reduces to 

wm 2 + vm 3 = —umi + wm 3 = —vm\ — wm 2 = 0. 

Since m 1? m 2 , m 3 > 0, the signs of u, v, and w must be such that 

(21) (i) u, w > and v < or (ii) u, w < and v > 0. 

Notice that from the first equation of each of the linear system in ( f20l) . 
as well as from the third equation of both systems, we can, respectively, 
conclude that 

§. <L 
m 2 = n 6 „ and m 2 ~ 



c 

b v b 



Then = Hz ^ H , a condition that is independent of 8 > 0. Two similar 
conditions follow by circular permutations. Therefore the systems in (120]) 
have the same set of solutions if 

, , a — c u — w b — a v — u c — b w — v 
22 — — = , = , = . 

b v c w a u 

We will further prove that these conditions are simultaneously satisfied 
only when a = b = c and u = v = w, a set of solutions which corresponds 
to equilateral triangles and, as a result, to equal masses. 

Let us first assume that the acute triangle is isosceles and not equi- 
lateral. Then, without loss of generality, we can choose a = c 7^ b and 
u = w 7^ v. The first equation in (T2"2l) is satisfied, whereas the other two 
take the same form, namely = which implies that av = bu. Since 
a, b > 0, it follows that u and v have the same sign, a conclusion which 
contradicts (T2Tj) . 

The last case to consider is that of the scalene acute triangles that are 
not isosceles. Without loss of generality, we assume that < a < b < c. 
Then 

a — c b — a c — b 

— — <0, >0, >0, 

b c a 

inequalities which, via ( 12"2"|) . imply that 

. . u — w v — u w — V 
(23) < 0, > 0, > 0. 

V w u 

But relations ([21]) provide only four possible ways to order u,v, and w, 
namely: (1) v < < u < w, (2) v < < w < u, (3) u < w < < v, and 
(4) w < u < < v. It is easy to see that, in each case, the order and the 
signs of these constants contradict at least one of the inequalities given 
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in f!23p . so scalene non-isosceles triangles cannot form nomographic orbits 
either. This remark completes the proof. □ 

7. Geodesic scalene triangles 

In this section we restrict our study to the case k > and consider 
a situation that was not captured in the previous section, namely the 
relative equilibria that rotate on great circles of a sphere generated from 
fixed points of the equations of motion. Such orbits do not exist for n < 
because, as shown in [7J and [8], there are no fixed points when the bodies 
move in hyperbolic space. 

Without loss of generality, we will analyze these orbits when the great 
circle is the equator, z = 0, and the rotation takes place around the z 
axis. This case is not captured by the system of differential equations 
given in equations (1151) and f[T7|) because one denominator cancels when 
the bodies reach the equator. When r is constant, i.e. the solution is just 
a relative equilibrium, instead of a homographic orbit with both rotation 
and expansion and/or contraction, the motion on the equator can be 
studied separately. In this case, no cancelling denominators show up 
when the motion takes place on the great circle z = 0. Nevertheless, no 
two bodies can be antipodal because the corresponding configuration is a 
singularity of system ([8]) (see [7J or [9] for more details), and therefore the 
motion doesn't exist. But acute triangles moving on the equator don't 
have antipodal bodies at their vertices, so singularities do not affect them. 

We will first focus on fixed points lying on the equator, and show why 
for every acute triangle there exist masses that provide a fixed point for 
the equations of motion if the bodies are placed at the vertices of the 
triangle. Then we explain how relative equilibria can be generated from 
fixed points. By providing a large class of counterexamples, we also prove 
that not any three masses generate fixed points, and consequently not any 
three masses can form relative equilibria that move along the equator. 

We can now state and prove the following result. 

Theorem 4. For any acute triangle inscribed in a great circle of the 
sphere S^, there exist bodies of masses m 1 ,m 2 ,m 3 > that can be placed 
at the vertices of the triangle such that they form a fixed point of system 
([8]) for n = 3 and n > 0, i.e. for the equations of motion of the curved 
3-body problem in the case of positive curvature. 

Proof. As mentioned earlier, we can assume that the great circle is the 
equator z = 0. To form a fixed point of system ([8]), the initial conditions 
of the bodies of masses mi,m 2 , and m 3 must satisfy, at the initial instant 
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t = 0, the constraints 



q*(0) = q*(0) = 0, z = l,2,3. 

Let qi(0) = (xi,yi,0),i = 1,2,3, be the initial position of the body of 
mass rrii, i = 1, 2, 3, on the equator z — 0. Using equations (|HJ), a straight- 
forward computation shows that the above conditions reduce to solving 
the linear homogeneous algebraic system 



(24) 



qi 2 m 2 + <?i 3 m 3 = 
9i 2 m 2 + 9i3™ 3 = 
q 21 m 1 + q 23 m 3 = 
q 21 m 1 + q 23 m 3 = 
qsimi + 932™2 = 

931™1 + 932™2 = 0, 

for mi, m 2 , and m 3 , where, for i, j — 1, 2, 3, i ^ j, 



KXiXj + KyiUj. 



Ql] (1-a?.) 3 / 2 ' % ' (1-a?-) 3 / 2 '^ 

But the first and the second equation of the system are linearly dependent. 
Indeed, multiplying the first equation by KX\, the second equation by Kyi, 
and adding the two new equations, we obtain an identity. Similarly we 
can prove the linear dependence of the third and fourth equation, as well 
as of the fifth and sixth equation. Therefore system ( |24|) can be reduced 
to the linear homogeneous system 



(25) 



qi 2 m 2 + g 13 m 3 = 
q 2 imi + g 2 3^3 = 
931^1 + q 32 m 2 = 0, 



in unknowns mi, m 2 , and m 3 . 

We will further show that this system has positive solutions. To achieve 
this goal we will first prove that detA = 0, where A is the matrix that 
defines system (125]) . namely 



Notice first that 



.4 



detA 





921 
931 



912 



932 



913 
923 





912923931 + 913921932- 

But qi2 and q 2 i have the same denominator, which is never zero. The 
same is true for the pair g 2 3 and g 32 , as well as for the pair g 31 and g 13 . 
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Therefore to prove that det A = 0, it is enough to compute the numerator 
E of det A. Notice that E has the form 

E = (x 2 -a 1 2Xi)(x 3 -a23X2)(xi-ai 3 x 3 ) + (x 3 -a 13 x 1 )(xi-ai2X2)(x2-a 23 x 3 ). 
A straightforward computation leads to 
E — 2(1 — 012013023)^1^2^3 + (012023 - ai 3 )xiX 2 + (013023 - ai2)x\x 3 
+ (ai 2 Oi3 - a 23 )x 1 xl + [a 13 a 12 - 023)^1X3 
+ (013023 - 012)^2X3 + (ai 2 o 2 3 - a i3 )x 2 xl. 

Without loss of generality, we can assume that x\ = 0, i.e. the body of 
mass mi is fixed at the point of coordinates (xi,yi, Z\) = (0, ft _1 / 2 ,0) on 
the sphere Sf. (see Figure [p. Then E = Sx 2 x 3 , where S = —013X3 + 
O13O23X2 + 012023X3 — ai 2 x 2 . Since x 2 , x 3 7^ 0, to prove that det A = it is 
enough to show that S = 0. Using again the fact that xi = 0, and writing 
the constants explicitly, we obtain 

S = -Kyiy 3 x 3 + n 2 y 1 y 3 (x 2 x 3 + ^2^3)^2 + n 2 ym(x 2 x 3 + y 2 y 3 )x 3 - Ky l y 2 x 2 

= -Kx 3 y x y 3 - Kx 2 y\V2 + Kx 3 y x y 3 (Kx\ + ny 2 2 ) + nx 2 y\y2{^\ + Ky^). 

But since kx\ + ny\ = kx\ + ny 2 = 1, it follows that 5 = 0, therefore 
det A = 0, so system (|25|) has other solutions than the trivial one. 

We must still prove that among these nontrivial solutions there is a class 
of strictly positive solutions. For this purpose, we will analyse the signs 
of the constants qij,i,j = 1, 2, 3, and see how they determine the signs of 
mi, m2, and m 3 . We already assumed that xi = 0. Since the triangle is 
acute, we can further suppose, without loss of generality, that x 2 > and 
X3 < 0, which means that the body of mass m 2 lies in the fourth quadrant 
on the equatorial circle, while m 3 is in the third quadrant (see Figured]). 
There are no other possibilities because the triangle is acute. (In fact we 
showed in [7] and [8] that bodies lying in the same hemisphere cannot 
form fixed points, thus excluding right or obtuse triangles as fixed-point 
candidates. Right triangles are also excluded because they form singular 
configurations.) 

Under the above assumptions, we will show that 

(26) <?i2,<?2i,?32 > and q 13 , q 31 , q 23 < 0. 

For this purpose, notice first that the denominators of g^, i, j = 1, 2, 3, are 
all positive, so we have to determine only the signs of their numerators. 
Also remark that the angles a± 2 and a± 3 from the centre of the circle 
corresponding to the arcs m\m,2 and mim 3 , respectively (see Figure [1]), 
both angles taken to be smaller than it, are larger than tt/2. Since 012 = 
cosai2 and a r3 = cosai 3 , it means that ai 2 , ai 3 < 0. Then, using the fact 




Figure 1 . The fixed point solutions formed by the masses 
rai,m2, and m 3 on the geodesic (equator) z = 0. 

that xi = 0,x 2 > 0,x 3 < 0, we can conclude that g i2 > and g i3 < 0. 
Since a±2 = a 2 i and a 3i = ai 3 , it also follows that g 2 i > and g 3i < 0. 

To prove the last two inequalities in fliZo 7 ]) . let /3 be the angle from the 
centre of the circle corresponding to the arc m 2 m 3 , a the similar angle 
corresponding to the arc m^A (see Figure d]), and 7 the similar angle 
corresponding to the arc m 2 A, all taken to be smaller than 7T. Then, 
obviously, a = /3 + 7. Also notice that x 3 = cosa,x 2 = cos 7, and 
a 23 = cos f3. So 

£3 — 0,23X2 = cos a — cos cos 7 = — sin (3 sin 7, 

which is negative because < 7 < 7r. Therefore g 23 < 0. Finally, 

£2 — ^23^3 = cos 7 — cos /3 cos a = sin a sin /5, 

which is positive because < a, (3 < ir. So g 23 > 0. To see now that 
system ( |25|) has positive solutions, it is enough to notice that the two 
constants showing up in each of its three equations have opposite signs. 
This remark completes the proof. □ 

A direct consequence of the above result is the possibility to generate 
relative equilibria from fixed points. This fact stems from the action pro- 
duced by elements of the rotation group SO(3) on a fixed point. From the 
mechanical point of view this means that we can obtain relative equilibria 
if we apply initial velocities, of equal speeds, tangentially to the geodesic, 
all oriented clockwise or all counterclockwise. This obvious remark to- 
gether with Theorem H] prove the following result. 

Theorem 5. Consider the curved 3-body problem on the sphere S"^, given 
by system flH]) forn = 3 and n > 0. Then for any acute triangle inscribed 
in a great circle of the sphere, there exist masses m!,m 2 ,m 3 > and 
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Figure 2. The initial positions of mi,m2, and 777,3, which 
form an isosceles triangle on the geodesic (equator) z — 0. 

initial velocities such that if the point particles are initially placed at the 
vertices of the triangle, the corresponding solution is a relative equilibrium 
that rotates on the great circle. 

The converse of Theorem H] is false, which means that not any three 
point particles of given masses can form a fixed point on a great circle 
of S 2 . We will provide a class of counterexamples in the case of isosceles 
triangles. In fact, the following result identifies all masses for which there 
are no isosceles triangles that lead to fixed points of the equations of 
motion of the curved 3-body problem for k > 0. 

Theorem 6. Consider the curved 3-body problem on the sphere S%, given 
by system (jSJ) for n = 3, k > 0, as well as the masses m\ = 7772 =: M > 
and 7773 =: 777 > 0. Assume that the initial conditions are such that the 
triangle having these masses at its vertices is acute and isosceles, with the 
equal masses corresponding to the base. Then, for M > 4m, no isosceles 
triangle can form a fixed point. 

Proof. Let qi(0) = (xi(0), yi(0), 0), i = 1,2,3, be the initial position of 
the masses with the symmetries of an isosceles triangle as in Figure [21 i.e. 
such that 

£i(0) = x, x 2 (0) = -x, x 3 (0) = 0, 

yi(0) = y 2 (0)=y, 2/3(0) = -^ 2 , 

z 1 (0) = z 2 (0)=z 3 (0) = 0, 

with < x, y < K -1 / 2 . The conditions that the coordinates form a fixed 
point are 

Qi(0) = q*(0) = 0, i = 1,2,3. 
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Using these conditions in system (jSJ), denoting x := x(0),y := y(0) and 
asking that x(0) = y(0) = 0, a straightforward computation leads us to 
the relations 

M — Anmy 2 .. M — Anmy 2 
AK l l 2 x 2 y ' ^ 4/c 1 / 2 x?/ 2 
Consequently 3 bodies lying at the vertices of an isosceles triangle form 
a fixed point of the equations of motion if and only if 

(27) M = 4Kmy 2 , with < y < k~ 1/2 . 

Then some necessary conditions that an isosceles triangle forms a fixed 
point are 

M 1 
0< < -. 

Anm k 

Since k > 0, the first inequality is always satisfied. The second inequality 
reduces to M < 4m. The corresponding x > is then obtained from 
kx 2 + ny 2 = 1 and leads to the same condition. So, for M > 4m, there 
are no isosceles triangles that form fixed points on a great circle of S 2 .. 
This remark completes the proof. □ 

Remark. When the triangle is equilateral, the evaluation of simp, where 
cp is the angle between the abscissa and the radius of the circle to rri2 (see 
Figure 

E}, shows that y = 1/(2k^ 2 ), so by we can conclude that 
M = m, in agreement with what we know about Lagrangian solutions. 
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